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We study a generalized branching random walk where particles 
breed at a rate which depends on the number of neighboring particles. 
Under general assumptions on the breeding rates we prove the exis- 
tence of a phase where the population survives without exploding. 
We construct a nontrivial invariant measure for this case. 

1. Introduction. Scientists have been studying models for the evolution 
of a population since the end of the 19th century, starting from the branch- 
ing process introduced by Galton and Watson in 1875 [5]. The need for 
more realistic models has led to the introduction of a spatial structure: the 
branching random walk and the contact process (briefly, BRW and CP, resp.) 
are perhaps the most natural generalizations. In the BRW model each indi- 
vidual has a fixed position on a connected graph, for example, the integer 
lattice "L d , and an exponential lifespan of parameter 1 during which it breeds 
on neighboring sites according to a Poisson process of intensity A > 0. The 
number of individuals allowed per site is unbounded. Requiring that a site 
can be occupied by at most one individual, one obtains the CP. Both these 
processes exhibit two possible behaviors: starting from a finite population, 
either the population faces almost sure extinction (subcritical behavior), or 
it survives with a positive probability (supercritical behavior). In the super- 
critical case the BRWs population grows indefinitely and the mean density 
of the population diverges. For the contact process, obviously there is no 
divergence of the mean density of the population because this quantity is 
a priori bounded. In the supercritical phase, the CP has two invariant ex- 
tremal measures (see [8]). It is known that there exists a critical value of 
A separating the two behaviors: if A is smaller than the critical parameter, 
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the process exhibits the subcritical behavior, while, for larger A, it exhibits 
the supercritical one. We denote by Abrw and Acp the critical parameters 
of the BRW and of the CP respectively. 

The observation of natural environments suggests to remove any a priori 
bound on the number of individuals allowed per site and to introduce a self- 
regulating mechanism on the birth rates, which should provide a surviving 
though nonexploding population. Indeed, some ecological systems seem to 
be in a sort of equilibrium where the density of a population neither tends 
to zero nor to infinity. One may argue that we could be observing a sub- 
critical or supercritical system during a too short time span, nevertheless, it 
seems natural to try to translate into mathematical terms the competition 
for resources (see, e.g., the discussion in [7]). Other authors have introduced 
models for self-regulating populations. For instance, in the case of a pop- 
ulation living on a continuous and homogeneous space, Bolker and Pacala 
[2] studied a process where the death rates depend on the local density cen- 
tered on the father. A slightly different model was considered in [4] where 
the reproduction rate depends on the local density centered on the father. 
The main technical tools are moment equations and stochastic differential 
equations respectively. A different approach is carried out in [3] where the 
population has no spatial structure and each individual can be affected by a 
gene mutation at birth; the evolution is studied as a Markov process in the 
trait space. 

We introduce a self-regulating mechanism where the birth rate is a de- 
creasing function of the local density at the location where the offspring 
would live. Moreover, noting that the spatial structure of the interaction 
between individuals in a biological population might be irregular, we study 
a population on a discrete (possibly nonhomogeneous) space. To this aim, 
we consider the following model, which we call restrained branching random 
walk (RBRW briefly). Consider an infinite connected graph X with bounded 
geometry (i.e., the number of neighbors of the vertices is bounded, e.g., 7j d ) 
as the environment where the population lives and let r](x) be the number 
of individuals living at the site x € X . The lifespan of each individual is 
an exponential random variable of mean 1. During its lifetime each indi- 
vidual tries to reproduce following a Poisson process of intensity A. Every 
time the clock associated to the Poisson process rings, the individual tries to 
send an offspring to a randomly chosen target neighboring site. The target 
neighboring site is chosen using the transition matrix P= (p(x,y)) x ^ y ^x of 
a nearest neighbor random walk on X, for example, the simple random walk 
on Z d . Call the target site y. The reproduction on y is effective only with 
probability c(r/(y))/A, where c : N — > M + is a nonincreasing and nonnegative 
function with c(0) = A. In this case the population living at y increases by 
one individual, otherwise nothing happens. 
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Observe that the process described above is a Markov process and includes 
the BRW and the CP as special cases (c = A and c = Al{ }, resp.). The formal 
construction of this process is carried out in Section 3, where the existence 
of a Markov process {r]t}t>o with state space Q C N x is proven. In general 
$7 is smaller than N x because we can only consider configurations rj such 
that rj(x) does not diverge too fast when x goes to infinity (see Section 2 
for more details). We prove that {r]t}t>o has different behaviors depending 
on c(0), c(+oo) := lmifc^ +00 c(k) and on the transition kernel p(x,y) (see 
Proposition 4.1 for the complete statement). 

Proposition 1.1. Let P n = (p( n \x,y)) x>ye x be the nth power of the 
transition matrix P and {r]t}t>o be the RBRW described above. Let us define 

p := limsup^^ ^/p( n )(x,y) and 9 := lim^oo ^sup x Y.yP^iv^) (notice 
that 9>p): 

(i) //c(0) < 1/p, then lim^ +00 E, v [r]t(x)] = for any finite n G fl, x G X; 

(ii) Ifc(0) > A C p, then lim^ +00 W[r] t (x)] > for any r] G O\{0}, x G X 
and P^limsupf^oo i] t (x) > 0) > 0; 

(hi) If c(+oo) > 1/p, then lim^ +00 ¥, v [r]t(x)] = +oo for any n G 0, \ {0}; 
(iv) If c(+oo) < 1/0, then limsup t _ >+00 E^fry^x)] < +oo uniformly for 
any bounded rj x G X . 

The critical parameters p and # (and Acp as well) depend only on P 
(hence, in the case of the simple random walk, on the geometry of the graph); 
in particular, Abrw = 1//° ( see i^} and [1]). We discuss further details in 
Section 2.1. The proof of Proposition 1.1 is quite simple and essentially based 
on coupling techniques with the CP and the BRW (and on explicit estimates 
on the moments of the BRW with immigration; see Lemma 3.3). 

Notice that, given a bounded initial state rj G SI \ {0}, if c(0) is sufficiently 
large and c(+oo) is sufficiently small, then by (ii) of Proposition 1.1, the 
population has a positive probability to survive indefinitely, while by (iv), 
almost surely, it does not explode. This is the ecological equilibrium phase 
we are looking for. It is quite natural to wonder if there is a stationary 
distribution for the population in this case. We prove that this is the case 
(see Theorem 4.3 for the complete statement). 

Theorem 1.2. Let {f/j}t>o be the RBRW described above and assume 
that c is such that c(0) > Acp and c(+oo) < 1/9. Then there exists a non- 
trivial probability measure p on {Q.,B{Vt)) which is invariant for {r]t}t>o- 

We construct this invariant measure as a limit of invariant measures of 
processes where the number of individuals per site is bounded. 
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Important examples are the RBRW on the d-dimensional lattice and on 
the homogeneous tree of degree n + 1 (both endowed with the transition 
matrix of the simple random walk). 

If X = Z, d and P is the simple random walk, we have that p = 9 = 1. 
So (i) of Proposition 1.1 implies that the population dies out when c(0) is 
smaller than the death rate 1, while (iii) states that the mean density of 
the population explodes when c(+oo) is larger than 1. Moreover, the system 
can reach the ecological equilibrium if c(+oo) < 1 and c(0) > Acp (i.e., the 
critical parameter of the CP on Z d ). In this case the stationary measure p 
given by Theorem 1.2 is translation invariant. 

If X = T n+ i (the homogeneous tree where the degree of each vertex is 
n + 1) and P is the simple random walk, we have p = 2y/n/(n + 1) < 1 = 9, 
whence, to ensure ecological equilibrium we require c(+oo) < 1 and c(0) > 
Acp (i.e., the critical parameter of the CP on T n +i), while, for the almost 
sure extinction, it is sufficient that c(0) < (n + l)/2y/n. The stationary mea- 
sure p given by Theorem 1.2 is translation invariant in this case as well. 
Considering different random walks on T n+ i leads to different values for 9 
(see Example 5.1). 

One may wonder how the two parameters p and 9 come to surface: the 
analysis of the two examples above, shows that on general graphs the behav- 
ior of interacting particle systems can be different than on For instance, 
it is known that on some fast growing graphs there is the so-called weak 
phase (see, e.g., [11] and [9], Part I, Chapter 4 for the CP on trees, [15] and 
[1] for the BRW on graphs and [12] for the BRW on Galton- Watson trees): 
the population can survive by drifting to infinity and leaving eventually any 
site. 

On Z, d , Abrw = 1 and, in the subcritical phase, starting from a bounded 
r\ with at least one individual per site, with probability one, there is no 
extinction and the expected number of individuals at a fixed site is a bounded 
function of the time t. 

On a general graph the subcritical phase of the BRW is further subdivided: 
Abrw = 1/Pj but only if A < 1/9 one can ensure that the expected number of 
individuals at a fixed site is a bounded function of the time t (starting from 
a bounded rj with at least one individual per site). Indeed, if A G (1/9, 1/p) 
and rj is as above, then there are examples where the expected number of 
individuals at a fixed site diverges as t goes to infinity (e.g., if P is the simple 
random walk on a homogeneous tree: see Example 5.1). 

We give here a brief outline of the paper. In Section 2 we give the defini- 
tions needed in the sequel and we introduce the generator of the process. The 
construction of the RBRW is carried out in Section 3. Since the state space 
of this process is not locally compact, the classical approach of Hille-Yosida 
cannot be used: we follow the ideas of [10]. Some of the results we prove 
are obtained via a coupling argument (see Proposition 3.5) with particular 
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BRWs (with immortal particles). Furthermore, we give explicit estimates of 
some moments of these processes (see Lemma 3.3). In Section 4 we prove 
our main results: Proposition 1.1 and Theorem 1.2 (see Proposition 4.1 and 
Theorem 4.3, resp.). Section 5 is devoted to final remarks, examples and 
open questions. 



2. Preliminaries. 



2.1. Graph geometry and random walks. Let X be a connected, nonori- 
ented graph, with bounded geometry (i.e., the number of neighbors of a 
vertex x, called degree of x, is uniformly bounded on X); denote by D the 
maximum degree of vertices on X. Let P = (p(x,y)) x , y ex be a stochastic 
matrix (although one can apply easily our methods to a substochastic P) 
such that p(x, y) > if and only if x and y are neighbors (we write x ~ y in 
this case). For any A C X, let A° := {x £ A : Vy ~ x, y G A} be the interior 
of A and let 

PA{x,y):-^ otherwise. 

The two parameters p and 9, associated to P, play a crucial role in dis- 
tinguishing between different behaviors of the RBRW (see Proposition 1.1). 



Recall that in Section 1 we defined p := limsup n ^ 00 yp^(x,y): this is usu- 
ally called the convergence parameter and it is independent of x,y G X (see 
[14])- 

Consider the space of (infinite) matrices endowed with the norm \\A\\ := 
su Px J2 y \ a xy\ ■ As usual, each matrix with a finite norm can be identified 
with a linear continuous operator from l°°(X) into itself. Let 9(A) be the 
spectral radius of the operator A; note that 9(A) = lim n ^ 00 

||4 n || 1/n (see, 

e.g., [13], Theorem 18.9). The parameter 9, defined in Section 1, satisfies 
9 = 9(P T ). The estimate of 9 is easy in some cases: 

(a) if there exists v. X i— ► [0, +oo), u ^= 0, such that u(x) < Cv(y) for any 
x, y G X (and some constant C > 0) and v(x) > J2 y v (v)p(V-, x ) [resp. v(x) < 
J2 y "(y)p(y, x)], then 9 > 1 (resp. 9 < 1); 

(b) if lim n ^ +0O (#{y : \y\ < n}) 1 /™ = 1, then 9 < 1. 

Note that in general p < 1 and p < 9. Using the above result (a), given 
a graph X, if P is strongly reversible (see [6]), then the amenability of X 
implies p = 9 = 1, while nonamenability implies p < 9 = 1. Roughly speaking, 
amenable graphs are graphs where the boundary of a finite set may be 
arbitrarily small if compared to the size of the set itself (see [16], page 112). 
Hence, p = 9 = 1 in the case of the simple random walk on Z rf , while there 
are examples where 9^1 and 9 > p, for instance, on homogeneous trees 
T n+ i [see Example 5.1 with p ^ l/(n + 1)]. 
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2.2. Configuration space. Following [10], fix a reference vertex xq G X 
and denote by |x| the graph distance between x and xq. Define a strictly 
positive function a : X — > R + by a(x) = M~\ x \ where M > (D — l) 2 . By this 
choice of M, for any z > 1/2, J2x a ( x ) z < +°° an d 

(2.1) £ g (x,y)a(i,)<Ma(a?) 

y 

for any substochastic matrix Q. Given r\ : X — > N, define ||t/|| := f?(a;)Q;(a;). 
The configuration space is f2 := {77 G such that 1 1 77 1 1 < +00}, while Qa '■= 
{r?:N A such that \\t)\\ < +00}. Note that the finite configurations, that is, 
the configurations ij Gil such that J2x v( x ) < +00, are dense in Q with this 
norm; moreover, the Borel c-algebra induced by the norm is the same as 
the one induced by the product topology. We introduce the usual partial 
order on £2, that is, £ < r\ if < r\(x) for any x £ X. We denote by 
the configuration identically equal to 0, by 1 the configuration identically 
equal to 1 and by 6 X the configuration which is equal to at any site but 
x, where it equals 1. We say that a function / : Q — > R is nondecreasing if 
£ < r] implies /(£) < f{rf). Given fi, v probability measures on £1, we say that 
v stochastically dominates fi and we write [i < v if for any nondecreasing 
function / we have fj,(f) < v{f) [where fj,(f) = J n fdfj]. 

From (2.1) we derive a useful bound on the transition kernel of the con- 
tinuous time random walk associated with P and with jump rate A > 0. 
Indeed, let 

n=0 

then the iteration of (2.1) gives 

(2.2) ^^(x,y)a(y)<e At ( M - 1 )a( a ;). 

y 

For A C X , we will also denote by 

(2-3) ^(^) : = e ~ A 'E^AW). 

77-. 

Clearly, p^ A (x,y) <p^(x,y), hence, the bound in (2.2) holds for these "re- 
stricted" kernels as well. 

2.3. Dynamics. Denote by Lip(O) the set of the Lipschitz functions on 
f2, and given / G Lip(O), let L{f) be its Lipschitz constant. For any / : — > K 
and xeX, define (d~f)(r]) = t[i ;+00 )(v(x))[f(ri-S x )-f(r])} and (d+f)(r)) = 
f(v+S x )~f(v)- Note that |(9±/)(r?)| < L(f)a(x) for any / G Lip(n), x £ X. 
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Fix a nonincreasing function c : N — > M + , a transition matrix P on X and 
define £ : Lip(Q) -> M Q by 



(2.4) (£/)(r?) :=£>(x) 



(0*/)fa) + E^fo)M*,i>)(3j7) fa) 



It easy to check that |(£/)(t/)| < L(/)[(c(0)M + 1)||??||], hence, £ is a well- 
defined operator on Lip(fi). 

3. Construction of the process. The main result of this section concerns 
the construction of a process having generator C given by (2.4). It is a stan- 
dard fact of the theory of countable state continuous time Markov chains 
that there exists a unique Markov process {rjt}t>o with generator given by 
(2.4) starting from any finite r/Gfi. The extension of this construction to 
more general configurations requires more sophisticated techniques. Our ef- 
forts in this direction may be summarized in the following proposition whose 
proof is the consequence of several intermediate steps. 

Proposition 3.1. There exists a unique semigroup {St}t>o °f opera- 
tors, St : Lip(fi) — > Lip(f2) such that: 

(i) {S t f)(v)=mf(Vt)] for f € Lip(fi) and r, finite. 

(ii) L(S t f) < L(/)exp(c(0)Mt) for f G Lip(fi). 

(iii) (S t f)( V ) = f( V ) + Jl{CS u f){ri) du for f £ Lip(fi) anrf^JI. 

(iv) Lei p. be a probability measure on O suc/i i/mi is finite. Then \x 
is invariant for {St}t>o [i.e., for any t > 0, fJ>[Stf] = fJ-[f] for any f £ Lip(f2)/ 
if and only if n[Cf\ = for any f 6 Lip(fi). 

Given Proposition 3.1, in order to define the process {r]t}t>o starting from 
any ?)£(!, according to [10] (see discussion after Theorem 1.4), one shows 
that {S*(}t>o can be extended to any measurable function / on U which 
satisfies either / > or \f(r))\ < C(l + \\r]\\) for some constant C. Thus, it 
identifies a unique Markov process on O (the RBRW) which we still denote 
by {vt}t>o- 

We start by constructing the process on a finite subset A C X. We need 
an auxiliary process defined on fi. Fix 7 > 0, c : N — ► M + , fcGN and define 
C?A : Lip(fi) -» R n by 



(5a/)(»7): = E 

(3.1) 



7 l A (x)(7?( a ;)-A ; ) + (9-/)(r ? ) 

+ ^(z) ^Pa(^, y)c(v{y))(dy f)(ri) 



8 



D. BERTACCHI, G. POSTA AND F. ZUCCA 



It is obvious that Q\ generates a Markov process {n t }t>o defined on 17a- 
Clearly, this process may be thought as a process on f2 where the particles 
outside A are "frozen" in the initial state. Furthermore, if this process starts 
from rjo G f2 such that r/o > kl, then obviously jjt > kl for any t > 0; in this 
case we say that there are k immortal particles per site. 

3.1. BRW with immortal particles. The estimate of the first and the 
second moments of the process generated by (3.1) will follow from a coupling 
with the BRW with k immortal particles, that is, the process where c(-) = A. 
Hence, we take a technical detour and study this particular process [or, 
equivalently, the BRW with constant immigration rate; see (3.4)]. 

Although in this section we are treating only finite sets, the following 
lemma is needed also in the countable case (see Remark 3.11). 

Lemma 3.2. Let Q be a (possibly infinite) matrix with q(x,y) G [0, 1] and 
finite norm and define qt(x,y) := e~ xt Y^=o{Xt) n q^ n \x,y) /n\. Let 
f : [0, +oo) — > l°°(X) be such that limt_ +00 f(t) = v. Given the system of linear 
differential equations, 

(3 2) = X (j2 ( i( x >y) u ( t iy) - u (t> x )) + (3u{t,x) +f{t,x) Vz, 

(u(0,-) = <p{-)? 

where (p G l°°(X) and [3 > A(l + @(Q), the corresponding solution satisfies 
lim t ^ 00 u(t, x) = ((A - [3)1 - \Qy x v for any x£X and ip G l°°(X). 

Proof. The proof is standard and we just sketch it. One can solve the 
system by considering the (stronger) Cauchy problem in l°°(X) 

(u(t) = -Au(t)+f(t), 
W0) = tp, 

where A = (A — /3)I — XQ. By our hypotheses we have that Re(cj(^4)) > 
e > 0, hence, \\e~ M \\ 0, || / * e~ As ds\\ < J +o ° \\e- As \\ds < +oo for any 
t G [0, +oo] and A is invertible. The solution is given by the well-known 
formula u(t) = e~ At (p + Jq e~ A ^~ s ^ f(s) ds. The first term tends to zero, while 
the second one can be written as J io e~ A ^~ s ^ f(s) ds + f tQ e~ A ^~ s ^ f{s) ds and 
the claim follows choosing to such that for t > to, \\f(t) — v\\oo is sufficiently 
small. □ 

Lemma 3.3. Let A C X be finite. Fix 7 > and k G N. Consider Q\ 
defined in (3.1) with c(-) = A > 0. Let {r]t}t>o be the process generated by 
Q\ starting from r\ bounded, r]>kl. Moreover, if A < 7/6, then there exists 
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two nonnegative constants Ui iA (k, A, 7) and U2 t \{k, A, 7) such that, for any 
x £ A, we have that 

(3.3) lim K>[Tk(x)] < U 1>A (k,X,l), lim E r '[( Vt (x)) 2 } < U 2>A (k,Xn), 
where the limits are attained uniformly with respect to x. 

Proof. Define {£t}t>o as & := rjt — kl. This process is a Markov process 
(viz., it is the branching random walk with constant immigration rate \k) 
and its generator is 



(3.4) (Wa0)(0:=£ 

x<=A 



y 

Obviously for any r/ G 0, such that rj>kl, we have 
E r >[r ]t (x)]=E r '~ kl Mx)]+k, 

(3.5) 

m(Vt(x)f] =E^ fel [(6(x)) 2 ] +2kE^ kl Ux)] + A: 2 . 

Choose £ G r2 and let m(t,x) := E^[^(x)], for any x G A; by basic semi- 
group properties we have that ^m(t,x) = E^Ha^X^)] (where 7r x is the 
projection on the x coordinate). By computing explicitly TL a -k x , we obtain 
that m satisfies the system (3.2) with Q = P T , f(t,x) = k\J2 y PA(y-,x) and 
[3 = A — 7. The claim follows from Lemma 3.2. 

To prove the second moment assertion, consider C(t, x, y) := E^[^ t {x)^t{y)], 
for any x,y G A. Using the same arguments as before, we obtain that C is 
the solution following system of linear differential equations: 

(3 6) i ^C{t,x,y) = -AC(t,x,y) + f(t,x,y) \fx,yeA, 
\c(0,x,y)=S(x)S(y), 

where 

A = 2( 7 I -XB), 
5((x,y),(xi,yi)) = d y {y 1 ) Vb x (xi) - , 

/(*, x, y) = Xk[ m(t, x) £p A (z, y) + m(t, y) ^p A (z, x) 



+ 5 x (y)X yk^p A (z, x) + ^2p A (z, x)m(t, z)j 
+ 5 x (y)>ym(t,x). 

The system (3.6) is formally equivalent to the one in (3.2) with X x X in 
the place of X. The results just obtained for m ensure that / satisfies the 
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assumptions in Lemma 3.2. Moreover, 9(B) < 0, since b^ n \(x,y), (x\,yi)) = 
2^" YJk=Q (k)P^ ( x i ' x)p ( - n ~ k ^ (yx , y) and the claim follows. □ 

Remark 3.4. It is known that, given the equations (3.2) with x £ A, 
where A is finite, an explicit expression of the solution, for any tp G f2, is 

u(t, x) = Qt(x, yMy) + £ f ^ {t ~ s) Qt-s (x, y)f(a, y) ds. 

y y Jo 

This formula represents the solution also for infinite A under mild assump- 
tions on /: suppose, for instance, that / > 0, J2 X f(t> x ) a ( x ) < +°° f° r some 
t > and dtf(t,x) is bounded on any compact set, uniformly with respect 
to x G A. Moreover, if we consider two families {Qa}a and {/a}a which are 
nondecreasing with respect to A and we denote by {ua}a the corresponding 
solutions, then we have that u\ f ux as A j X [hence, one may replace the 
upper bounds in (3.3) with U\ x and Uzx which clearly are uniform with 
respect to A]. 

We note that the bound (3.3) can be obtained starting from any 77 £ CI (not 
necessarily bounded) by using similar computations, if A < 7/||-P r || < j/Q- 
It is easy to show that \\P T \\ = 1, for instance, for any symmetric random 
walk. 

3.2. The finite volume process. The following proposition shows how to 
construct a monotone coupling of different processes generated by (3.1). 

PROPOSITION 3.5. Fix N G N and let Ai C • • • C An C X be finite sub- 
sets. Fix ki, ■ . ■ , G N ; 71, . . . ,7jv G [0, +oo) and let c\, . . . ,cjv : N — ► [0, +oo) 
6e nonincreasing functions. Then for any fixed (770,1 , • • • , 770 jv) £ stic/i i/iai 
> ^ftl /or fa 6 {1, ••• ,AT}, i/iere exists a Markov process 

{(?7t,i, . . . , ??i,/v)}t>o on swc/i i/iat, /or any h £ {1, . . . jiV}, i/ie semigroup 
associated with the process {r]t,h}t>o has generator Q\ h . Furthermore, as- 
sume that k h < k h+1 , j h > j h+1 , c h (k h+1 + n) < c h+ i(k h+1 + n) for any 
n G N and r] 0th < ??o,h+i for any h £ {1, . . . , N - 1}. Then n t .\ < ■ < rjt,N 
for any t>0. 

Proof. It is enough to consider the processes on CIa n . Choose (770,1, . . . , 
77o,/v) m (CIa n ) N such that 770^ > k^l for any h G {1,...,N} as the ini- 
tial configurations. For any x G An, let A(x) := max{ 770,1 (x), ... ,rjo t Ar(x)}, 
7 := max{7i, . . . ,7tv}, c := max{ci(0), . . . , cjv(0)}. Choose an independent 
family of exponential clocks, two per site x G A at: one of parameter jA(x) 
which controls the deaths and one of parameter cA(x) which controls births. 
Define (77^1, . . . , r]t,N) '■= (Vo,i, ■ ■ ■ > Vo,n) for any t < r, where r is the time of 
the first ring of the collection of clocks. Assume that the clock which rings 
first is at site x: 
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• If the clock is a death clock, then for any z 7^ x put rj Tt h(z) := r] T - t h(z), 
pick a uniform U in the interval (0, 1) and define r) T h{x) := 7? r _ h(x) — 1 for 
any h such that U < (jh(r]r~ ,h( x ) ~ k h) + )hA(x) and r) T>h (x) :=r] T - ih (x) 
otherwise. Finally, restart the procedure from (r/ Ti i, . . . ,t] t ,n)- 

• If the clock is a birth clock, then for any z 9^ x put t] T h(z) := Vr-,h( z )- 
Choose at random, accordingly to the transition matrix P, a site y among 
the neighbors of x. Now pick a uniform V in (0,1) and define r) T h(y) '■= 
i] T - t h(y) + 1 for any given h such that V is not larger than r] T ~ h (x) 
x Vh h {x,y)c h {r] T - )h {y))/{cA(x)p(x,y)) and T] Tjh (y) = T] T -, h (y) otherwise. 
Finally, restart the procedure from (777-1, . . . ,rj T ^). 

It is a simple exercise to check that this construction leads to the desired 
coupling. □ 

In the remaining part of this section we prove some basic bounds on the 
semigroup {St t \}t>o generated by Q\. We need these bounds to extend the 
construction of the process to an infinite A C X. The next result shows that 
the semigroup {St t \}t>o generated by Q\ maps Lip(O) into itself. The proof 
follows closely the proof of Lemma 2.1 in [10]. 

Lemma 3.6. Let Ac X be finite and {St t \}t>o be the semigroup gener- 
ated by Qa- Then for any f G Lip(f2), 

L(S t , A f)<L(f)exp{c(0)Mt). 

Proof. Take £ Q (£,£ > kl) and consider the monotone coupling 
{{vhvhVt,Vt)}t>o of Proposition 3.5 such that (t)q, rfi, rft, rfi) := (£A<,£,C, 
£ V Q. This means that ryj <r)t,Vt — Vt f° r an y t>0. Therefore, 

\(s tA f)(o - (su/)(oi = \nf{ri 2 t ) - fm\ 

<nf{ni)-f{4)\]<mnH-4\\\- 

To bound this last term notice that by monotonicity 

II??? - Vt\\ <J2 a ^ x ">^t( x ) - vlix)). 

X 

Furthermore, for any x 6 A, we claim that 

(3.7) J^t(x) ~ Vl(x)] < c(0)^p A (2/,x)E[ % 4 (y) - r, l t {y% 

y 

which implies, by (2.1), that 



-EOl^-^IO^c^MEtll^-^ll]. 
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This gives 

E\\\r£-r£\\]<E[\\rjb-T&\\]ewmm 

and the proof is complete. To obtain (3.7), use the generator of the coupled 
process or better notice (see the proof of Proposition 3.5) that the rate of 
the transition r)f(x) — T)j(x) — ► rjf(x) — T)\(x) + 1 is 

c ( r it- i x )) pa(j/j (y) ~ Vt— (?/))• 

y 

The following result is a simple consequence of Lemma 3.6. The proof is 
the same as the one of Corollary 2.5 in [10], hence, we omit it. 

Corollary 3.7. Let A C X be finite and {<St,A}t>o be the semigroup 
generated by Ga- For any f : Q — ► R such that \f(ri)\ < CjHr/H for all rj > kl 
and for some constant Cf > 0, we have that 

\(S t ,Af)(v)\<Cf\\ri\\exv(c(0)Mt). 

The next two results are the analogs of Lemma 2.6 and Lemma 2.7 in 
[10]. 

Lemma 3.8. Let AcA'cI be finite subsets. Fix 7 > and k G N. Let 
Ga,Ga' be defined by (3.1). Then for any f G Lip(O), n > kl, 

\(GA'f)(v)-(GAf)(v)\<L(f)^ + c(0)M)Y,U'\Ao(x)v(x)a(x). 

X 

Proof. The proof can be obtained by direct computation. □ 

Lemma 3.9. Let A C A' C X be two finite subsets. Fix 7 > and k G 
N. Consider the semigroups {St,A}t>o an d {St t A'}t>o associated with the 
generators Qa-iQa 1 defined by (3.1). Then for any f G Lip(f2) and n > kl, 

\(S t! A'f)(v)-(S t , A f)(v)\ 

< L(/)( 7 + c(0)M)e c (°) Mt 5:a(x)l A , VA o(x) ?? ( y ) f p«$(y,x) du. 

x,y Jo 

PROOF. Note that 

(S t ,A'f)(r]) - (S t ,Af)(v) = [ {S u ,a>(Ga> ~ GA)S t - u ,Af){r]) du. 

Jo 

By Lemma 3.8, 

\((0A> ~ GA)S t -u,Af)(v)\ < L(^- w ,a/)(7 + c(O)M) tA'\Ao(x)v(x)a(x). 
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By Lemma 3.6, 

L(S t ^ A f)<L(f)e c ^ M ^ u l 
Using this last estimate and the positivity of S u> a>, we get 
\(S u> a'(Qa' ~ G\)St-u,Af)('n)\ 

(3.8) 

<L(/)e c (°) M (*-")( 7 + c(0)M)^a(x)l A ,\ A o(x)5' w , A K7r a; )( 77 ). 

X 

By Proposition 3.5, (S Ut Anr x )(v) ^ E A ,[r/ u (x)], where {r/ t }t>o is the process 
generated by (3.1) with 7 = 0, c(-) = c(0) and k = 0. By Remark 3.4, we 
know that 

Ei / [ Vu (x)] = e^Y,p c u%(y^)v(y)- 

y 

Plugging this bound in (3.8), we get 

\(S u ,A'(Ga' - GA)S t -u,Af){r])\ 

< L (/)e c (°) M *( 7 + c(O)M) a(x)l A ,\ A o(x)^ A ) ,(y, x)r,(y), 

x,y 

which concludes the proof. □ 

3.3. Finite volume approximation. Following [10], we construct the pro- 
cess on X as a limit of processes defined on A finite. For any n G N, define 
A n := B(xo,n), that is, the ball of radius n and center xq. 

Proposition 3.10. Fix 7 > and k G N. For am/ n G N, consider the 
semigroups {StA„}t>o generated by Qa„ defined in (3.1). For any fixed t > 0, 
/ G Lip(fi) and 77 G fi, ry > kl, the sequence {St t A n :n £ N} is a Cauchy 
sequence. 

Proof. Assume that m < n, then by Lemma 3.9, 
\(S t ,A n f)(v)-(St,A m f)(v)\ 

< L(/)( 7 + c(0)M)e c ^ Mt £ a^l^c, (x)r/(y) /* p^(y, x) du. 

x,y J® 

We have to show that 

r c(0) 

can be dominated uniformly in n G N by a function cj>(x,u) G L X (X x [0,i]), 
where the measure on X is «(■)• The result follows by dominated convergence 



^p C u?A n (y> x )v(y) 
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since lim min ^. +00 Ia^ao, ( x ) — hm m _> +00 1(ao )c(x) = 0. We claim that we 

can take <j)(x,u) = T, y v(y)Pu (y,x), indeed, p c u { ^ n {y,x) < p$°\y,x) and by 
(2.2), we have 

a; y y 

= e c ^ M - 1 ^\\ V \\eL\[0,t\). □ 

The proposition above allows us to define for any t > 0, / £ Lip(r2) and 
//£!!, 77 > fcl: 

(S t /)fa):= lim (S ttAn f)(r}). 



Remark 3.11. It easy to show that with this definition we can drop 
the hypothesis that A is finite (take A f X) in Proposition 3.5, Lemma 3.6, 
Corollary 3.7, Lemma 3.8 and Lemma 3.9. The same can be done in Lemma 
3.3, since one proves that ¥, v [r] t (x)] = lim^x m(t, x) + k and K v [n t (x)r]t(y)] = 
limA-fx C(t,x,y), where the limit functions satisfy the corresponding differ- 
ential systems. Note that, in particular, the latter is not obvious, because 
Tj i — > r](x)r](y) ^ Lip(f2). Moreover, the process generated by (2.4) is mono- 
tone as a consequence of Proposition 3.5. 

Proposition 3.12. For any t>0, f £ Lip(fi) and rjefl, f]> kl, define 
(S t f)(n):= lim (S t , A J)( V ): 



1. {Sf}t>o is a semigroup. 

2. For all f £ Lip(fi), j] £ SI, 

(S t f)(r,) = m+ [\gS u f)(r,)du. 
Jo 

Proof. These properties can be proven exactly as in [10] (page 451 and 
Lemma 2.12) by using Lemma 3.6, Corollary 3.7, Lemma 3.3, Lemma 3.8 and 
Lemma 3.9 instead of Lemma 2.1, Corollary 2.5, Lemma 2.6 and Lemma 2.7, 
respectively. □ 

Among the properties of the semigroup {St}t>o which can be proven, we 
state the one which we need in the next section. 

Proposition 3.13. Let [i be a probability measure on £1 such that f^[\\rj\\] 
is finite. Then fi is invariant for {St}t>o (i-e., for any t>0, fJ>[Stf] = fJ-[f] 
for any f £ Lip(Q) ) if and only if fJ>[Qf] = for any f £ Lip(f2). 
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Proof. See the proof of Corollary 2.17 in [10]. □ 

Proof of Proposition 3.1. It is easy to show that (i) holds. The 
claim (ii) follows from Lemma 3.6 and Remark 3.11, while Propositions 3.12 
and 3.13 imply (hi) and (iv), respectively. □ 

4. Ecological equilibrium and invariant measure. In this section we study 
the behavior of the RBRW constructed in Section 3. In particular, Proposi- 
tion 1.1 and Theorem 1.2 are proven (see Proposition 4.1 and Theorem 4.3 
below). The main tool is the coupling between this monotone process and 
suitable contact and BRW processes. 

PROPOSITION 4.1. Let {r] t } t >o the RBRW generated by (2.4): 

(i) If c(0) < 1/p [resp. c(0) < 1/p], then limt^+oo rjt(x) =0 a.s. (resp. 
lim^+oo E^r/^x)] = 0) for any finite i] £ £1, x E X ; 

(ii) i/c(0) > A C p, then lim t ^ +00 E%t(a;)] > for any r] G \ {0}, x £ X 
and P^limsup^^ rft(x) > 0) > 0; 

(hi) 7/c(+oo) > 1/p, then lim^+oo E^f^x)] = +oo for any r] £fl\ {0}; 
(iv) If c(+oo) < 1/9, then limsup t ^ +00 E* 7 [?/i(x)] < +oo uniformly for 
any bounded ijell, x £ X . 

Proof. Recall that Proposition 3.5 and Lemma 3.3 hold for A = X (see 
Remark 3.11): 

(i) By Proposition 3.5, we can couple the process with a branching ran- 
dom walk {Ct}t>o starting from n with birth rate c(0) such that rjt < Ct- The 
first part of the claim follows by noting that Ct dies out almost surely (see 
[1], Theorem 3.1). As for the second part, the assertion is a consequence of 
Lemma 3.3 and Remark 3.11 [since k = 0, one can choose Ui t x(0, A, 1) =0]. 

(ii) By Proposition 3.5, there exists a supercritical site-breeding CP 
{Ct}t>o starting from n A 1 with birth rate c(0) > Acp and such that (t<rit. 
Theorem 4.8 of Chapter VI in [8] yields to the conclusion. 

(hi) By Proposition 3.5, we can couple {rjt}t>o with a branching random 
walk {Ct}t>o starting from r] with birth rate c(0) > 1 such that Q < and 
the claim follows. 

(iv) In this case there exists k € N such that c(k) < 1/9. By Proposi- 
tion 3.5 and Remark 3.11, there exists a process {(t}t>o generated by (3.1) 
with k = k, 7 = 1, birth rate c(k), such that Q > ijt. By Lemma 3.3 and 
Remark 3.11, we have 

limsupE^ 1 ^^)] < lim ¥7> vkl [( t (x)] < Ui, x (k,c(k),l). 
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□ 



Remark 4.2. The condition ]hnswp t ^ +00 W , [r]t(x)] < +00 implies that 
P r '(lim t _ t+00 r/ t (x) = +00) = 0, but r/ t (x), as a function of t, could be un- 
bounded almost surely. 

The remaining part of this section is devoted to the proof of the following 
theorem. 



Theorem 4.3. Let {St}t>o be the semigroup generated by (2.4) using 
Proposition 3.1. Assume that c(0) > Acp and c(+oo) < 1/8. Then there ex- 
ists a nontrivial probability measure /x on (f2,£>(f2)) such that fj,\Stf\ = n[f] 
for all t>0 and f E Lip(Jl). 



We need some preparatory results. 



Lemma 4.4. Assume that c:N^ [0, +00) is a nonincreasing function 
such that c(0) > Acp, while c(+oo) < 1/8. For any n G N, define c n := 
cl [0,71-1] an d consider the generator 



(C n f)( V ) :=5>(s 



(0-f)(v) + Y, c rMy))p(x,y)(d+f)(7 ] ) 



y 

Then there exists [i n probability measure on £1 such that: 

(i) fi n C n = 0; 

(ii) the sequence {fi n : n > 1} is nondecreasing with respect to the stochas- 
tic ordering of measures; 

(hi) denote by v\ the nontrivial invariant probability measure of the CP 
on X with parameter A := c(0) > Acp (see [8], page 265). Then v\ < \i n for 
any n > 2; 

(iv) the sequence {// n }neN is tight. 

Proof. Notice that C n is of the form (2.4) so it generates a Markov 
process {r]t <n }t>o: 

(i) The process {r] t . n }t>o is monotone because of Proposition 3.5 and 
Remark 3.11. If the initial condition is nl, then by standard arguments 
(see [8], Chapter III, Theorem 2.3), r/ tj „ =>• fi n as t — > +00. Furthermore, 

(ii) For any n > 2, by Proposition 3.5 and Remark 3.11, there exists a 
monotone coupling between {r]t t n}t>o, starting from nl, and {r]t, n +i}t>o, 
starting from (n + 1)1, such that rjt in < f]t n+i f° r an Y t>0. Let / : f2 — > R 
be a nondecreasing function, then E nl [/(7ft iri )] < E^" 1 " 1 ) 1 [f (i]t n+i)] f° r an Y 
t > 0. By taking the limit, as t — > +cxo, we get fJ> n (f) < / u n+i(/)- 
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(iii) By Proposition 3.5 we can couple {r]t.2}t>o, starting from 21, and a 
supercritical CP {£j}t>o, starting from 1, with parameter A = c(0) in such 
a way that & < 77^2 for any t > 0. 

(iv) Note that c(+oo) < 1/6* implies fc := inf{/c G N : c(k) < 1/9} < +oo. 
Take n >k and observe that by Proposition 3.5 and Remark 3.11 there 
exists a monotone coupling between {rjt t n}t>o, and the BRW {(t}t>o gen- 
erated by (3.1), with k:=k, 7 = 1, and parameter c(k), both starting from 
nl. Since r/ 4 , n < Ct for any t > 0, then E nl [rj t , n {0)} < E nl [Ct(0)]. By taking 
the limit as t —> +oo and using Lemma 3.3 and Remark 3.11, we have that 
H n (r](x)) < Ui t x(k,c(k), 1) for any n > k,x G X. Hence, there exists a con- 
stant C := Ui t x(k,c(k), 1) such that, for any r > and n > k, we have, by 
the Chebyshev inequality, fJ.n(v £ : ^(^O > r) <C/r. Let us fix ^4 > and 
define r(x) := A/y / a(x) for any xGl. We have 

/x n (77 G O : 77(3;) < r(x) for any x G X) 

= 1 — /i n (rj G : there exists i£l such that rj(x) > r(x)) 

> 1 - ^Vniv £ : Vi x ) > r ( x )) 

X 

^r(x) 

By our assumptions, J2 X V a ( x ) < whence, for any e > 0, we can choose 
A so that /j„(?j G : 77(2;) < r(x) for any x G X) > 1 — e for any n>k. The 
subset iT := {77 G 'v( x ) — r i x ) fo r an y x £ -^l °f ^ * s compact. In fact, 
K = ria:[0,r(x)i since 77 G Hr([p! r ( x )\ ^ implies that 

^7/(x)a(x) <A^2^a(x) < +00, 

that is, 77 G fi. □ 

Since the sequence {/J n }neN is tight and monotone and since the set of 
continuous, monotone functions separates the set of probability measures, 
then the sequence converges weakly to a probability measure on $7, say, /x. 

Moreover, /x inherits all the symmetries of X and P: if T is a bijection 
of X onto itself such that x ~ y if and only if Tx ~ Ty and p(Tx,Ty) = 
p(x,y), then /i(T~ 1 (A)) = fi(A) for all measurable set of configurations A. 
In particular, if P is the simple random walk on Z rf or on the homogeneous 
tree of degree n, then /x is translation invariant. 

By the previous lemma, /x n < /x for any n G N. Furthermore, < v\(r){x)) < 
fJ>(v( x )) — Ui t x(k,c(k), 1) (by the same bound on /x n ), hence, /i is not 5o and 
/4IMI] < +°°- We prove that /x is invariant by showing that fJ,[£f) = for any 
/ G Lip(f2) (see Proposition 3.13). In order to see this, we need a preparatory 
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lemma, indeed, in the proof of Proposition 4.6 we need that fi n (Cf) — > fJ-(Cf) 
as n — ► +00, but this does not follow directly from \i n =>■ [i because Cf is 
unbounded. 

Lemma 4.5. Let {/%}neN be a nondecreasing sequence of probability 
measures on Q and assume that \i n => fi as n — ► +00 and /i[|M|] < +00. 
For any m G N, 77 G $7 and g:^ — > R, define the configuration rj m (-) : = 
ls(i .m)(') J 7(') andg m {rj) :=g{rf n ). Assume that g:£l— >R satisfies the fol- 
lowing: 

(1) f/iere exisis C > suc/t i/iai | ^7(77) | < C(||n|| + 1) /or any 77 G fi; 

(2) /inllg — £f m |] — > asm-» +00 uniformly in n G N. 

T/ien fj, n (g) — > n(g) as +00. 
Proof. We have that 

IMs] - mMI < iMnb - 5 m ]l + MsH - + Hg m -g]\. 

By hypothesis (2) and the dominated convergence theorem, the first and 
last term on the right-hand side of the above inequality may be made small 
uniformly in n by taking m sufficiently large. Fix m G N such that these 
terms are smaller than e > for any n G N. For the middle term, define 
g m ' k (v) ■= 5 m (^)l(-oo,fc](b m («)l), k G N and write 

(4.1) 

< \vn[g m - g m ' k ]\ + \vn[g m ' k ] - n\g m ' k ]\ + Hg m ' k -g m ]\- 

Note that by hypothesis (1) and elementary bounds, 

<C(||r;||+l)i (fc/ cr-i I +oo)(hll) 

By monotonicity, the first and the last term on the right-hand side of (4.1) 
can be bounded above by //[t^]. Furthermore, hm/ c ^_|_ 00 v k = and v k (r]) < 
C(\\r/\\ + 1), so, by dominated convergence, the first and the last term on 
the right-hand side of (4.1) can be made smaller than e by taking k large. 
Finally, fix k G N large enough, and observe that the middle term on the 
right-hand side of (4.1) goes to as n — ► +00 by weak convergence. □ 

Proposition 4.6. Let /j be the weak limit of the sequence {// n }neN de- 
fined in Lemma 4-4> then n\Cf\ = for any f G Lip(f2). 

Proof. We start splitting 



(4.2) \fi[Cf]\ < M£/]-/i n [£/]| + Mn [|£/-£ n /|]. 
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Roughly speaking, the first one of these two terms goes to by weak con- 
vergence, while the second one goes to since C n f — > Cf. 

By Lemma 4.5, — H-n[£f}\ — > if we can show that condition (2) is 

satisfied by g := Cf [condition (1) is easily verified]. Observe that 

\(£f)(v)-(cf)(r)\ 

\(d-f)(v)-(d-f)(fT)\ 



x 

(4.3) 



+ c(o)J2p(x,y)\(d+f)(v) - (9+fW 



\(d-f)( v )\+c(o)J2p^y)Hf)W\ 



It easy to show that L(df f) < 2L(f) for any / G Lip(Q), x £ X. Thus, the 
first of the two sums on the right-hand side of (4.3) is dominated by 



2L(/)(1 + C(0))\\V ~ fTW E *-B(xo,m)(xW> 



Moreover, 



h-V m \\Yl 1 B(x ,m)(xMx) = Y^oi(y)t B ^ 0:m) (x)t B{xotm)C (y)rj(x)rj(y). 

x x,y 

This implies that the fj, n mean of the first term on the right-hand side of 
(4.3) is smaller or equal to 

2L(/)(1 + C(0)) E ^B(x ,m)(x)Z B (x ,m)z(y) a (y)Vn[v(x)v(y)}- 
x,y 

This term goes to uniformly in n G N since ^n[i]( x )'n(y)] — ^[vi^viv)] — C 
by Remark 3.11 and Lemma 3.3 while, by our choice of M > (D — l) 2 , 

E 1 B{x ,m) (x)^B(x ,m)V (v) a {v) '™ °" 

The second term on the right-hand side of (4.3) can be dominated by 
L(/)(l + c(O)M) Y, Z Bixo , m) c(xHx)a(x), 

x 

hence, its /i n mean converges uniformly in n to 0, since it is not larger than 

U hX (k, C(k),9) £ ^B(x ,raf (*)«(*) ^ 
x 

(again use Lemma 3.3). 
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We are left with the proof that the second term on the right-hand side 
of (4.2) goes to as n — > +oo. Observe that since p n is concentrated on 
{r/ : r] < nl}, 

p n [\Cf-C n f\] =c(n)J2p{ x iy)Vn[v(x)l {n} (r)(y))\d+f(r))\] 

< c(0)L(/) y)a(y)p n [r](x)t[n, +oo)(rj(y))}. 

By the Schwarz and Chebyshev inequalities, Lemma 3.3 and Remark 3.11, 
we have that 

[n,+oo) 

i c 

< sj p[i](x) 2 ]p(r](y) >n)<-j=^0, 

where C does not depend on x and y. □ 

Proof of Theorem 4.3. It follows from Lemma 4.4, Lemma 4.5 and 
Proposition 4.6. □ 

5. Final remarks and examples. In this paper we consider mainly local 
survival (i.e., persistence of the population in a fixed site). We already ob- 
served that, for general P (e.g., for the simple random walk on a general 
graph), a weaker type of survival (the weak phase) is possible for both the 
CP and the BRW. One can associate to this ylobal phase two critical pa- 
rameters (which coincide with Acp and AbrWj f° r instance, on Z rf ). Clearly 
this phenomenon could be observed also in the evolution of a RBRW. 

Proposition 4.1 does not describe the behavior of all possible RBRW: in 
particular, it is not clear what happens if 1/p < c(0) < Acp or if c(+oo) G 
(1/9,1/ p). As for this last question, one would ask whether this interval 
may be nonempty. We already noticed that p = 9 on amenable graphs with 
a strongly reversible random walk (such as the simple random walk on Z rf ), 
nevertheless, there are examples of graphs (see Example 5.1) where p ^ 9 
and for any A € (1/9,1/ p), the BRW starting from a finite configuration 
vanishes locally with probability 1, while starting from a bounded config- 
uration greater than 1, the expected number of individuals at a fixed site 
diverges. Roughly speaking, this is possible on graphs where the contribu- 
tion of far distant individuals is not negligible; indeed, if P is symmetric, 
this behavior is equivalent to the existence of a weak phase for the BRW. 
One may conjecture that the critical parameter of this phenomenon may 

be 1/limsupn^oo \Jj2xP^ n H x i v) (which does not depend on y). Finally, an- 
other open question is the extremality of the invariant measure p. 
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Example 5.1. Let us consider the homogeneous tree T n+ i where the 
degree of each vertex is n + 1 and choose n > 2. Fix a reference vertex o 
and p 6 [0, 1/n]. Given two neighbors x and y, we define p(x,y) as p if 
\x\ + 1 = 1 2/ 1 > 2 (recall that \x\ is the distance from x to o), as \ j{n + 1) if 
y ~ x = o and as 1 — np otherwise. By using standard generating function 
techniques and the fact that ||P T || > 9 > limsup,^^ \J^Z y p^ n \y,x) > p (for 
all x), it is easy to show that Note that p < 1 for all p > l/2n and that p < 9, 



Range of p 


P 









[0,l/2n] 


1 




1< n - (n 2 - l)p < 6 < (n + 1) ( 1 - 


np) 


(l/2n,l/(n + l)) 




-np) 


1< n - (n 2 - l)p < 6 < (n + 1)(1 - 


np) 




2vV(n 


+ 1) 


6 = 1 




[l/(n + l),l/n] 


2^/np(l 


-np) 


6 = n- (n 2 - l)p< 1 




for instance, if p = 


1/(71 + 1), 


that is, 


the simple random walk. In this 


case 



= 1; by using the explicit solution given in Remark 3.4 (with (p = 1), and 
noting that, by translation invariance, YjxP^ n \ x ^v) = \\(P T ) n \\ >9 n = l, we 
have that lim^oo E^f^ (x)] = +oo if A > 1. Moreover, considering the BRW 
with at least one immortal particle per site, if A > 1, then the expected 
number of individuals on a fixed site diverges as t goes to infinity. 
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